Superradiance mediated by Graphene Surface Plasmons 
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We demonstrate that the interaction between two emitters can be controlled by means of the 
efficient excitation of surface plasmon modes in graphene. We consider graphene surface plasmons 
supported by either two-dimensional graphene sheets or one-dimensional graphene ribbons, showing 
in both cases that the coupling between the emitters can be strongly enhanced or suppressed. The 
super- and subradiant regimes are investigated in the reflection and transmission configurations. 
Importantly, the length scale of the coupling between emitters, which in vacuum is fixed by the free 
space wavelength, is now determined by the wavelength of the graphene surface plasmons that can 
be extremely short and be tuned at will via a gate voltage. 



I. INTRODUCTION 

In the past few years, plasmonics has emerged as a way 
to control light at the subwavelength scale [THS]. Its po- 
tentiality is based on the properties of surface plasmons 
(SPs), which are surface electromagnetic (EM) waves 
coupled to charge carriers. These surface waves propa- 
gate along a metal-dielectric interface and are character- 
ized by their subwavelength light confinement and long 
propagation lengths. Recently, it has been shown that 
graphene 0], which has remarkable optical [5] and opto- 
electronic [5] properties, also supports the propagation 
of SPs. These EM modes bounded to a graphene sheet 
have been studied theoretically [7HT2] and very recently 
observed in experiments [T31 H3] • 

The properties of SPs in graphene have attracted great 
attention as they appear as an alternative for many of 
the functionalities provided by noble-metal SPs with the 
advantage of being tunable by means of a gate poten- 
tial. For instance, by designing spatially inhomogeneous 
conductivity patters in a graphene sheet one can have 
a platform for transformation optics and THz metama- 
terials 15 . Additionally, strong light-matter interaction 
between SPs and quantum emitters in graphene has been 
proposed based on the high decay rates of emitters close 
to graphene sheets [TH]. The field patterns excited by 
a nanoemitter in graphene were analyzed in Ref. 1171 
demonstrating high field enhancements and long prop- 
agation lengths for the SPs. Furthermore, fluorescence 
quenching in graphene has been proposed as a probe of 
the evidence of plasmons [TSJ QI5] . 

The paper is organized as follows: in Sec. II we in- 
troduce the theoretical formalism and characterize the 
coupling between an emitter and the GSPs supported by 
a 2D graphene sheet. In Sec. Ill we study the inter- 
action between two emitters mediated by the 2D GSPs. 
The presence of GSPs leads to interesting phenomena 
such as the super- and subradiant regimes, where the 



radiation properties of several emitters are enhanced or 
suppressed. We introduce the total normalized decay fac- 
tor, 7, which characterizes the superradiant regime, and 
discuss its tunability. We analyze the evolution of 7 with 
the vertical distance between the emitters and the sheet 
in the transmission configuration, and with the in-plane 
distance between the emitters in the reflection configura- 
tion. In Sec. IV we consider the coupling of two emitters 
by GSP in one-dimensional (ID) graphene ribbons |21j . 
We show that confining in ID leads to longer interac- 
tion ranges between the emitters. Moreover, we study 
how the coupling between the emitters is affected by the 
presence of the ribbon's edges. Finally, our main results 
are summarized in Sec. V. 



II. COUPLING BETWEEN AN EMITTER AND 
GRAPHENE SURFACE PLASMONS 

First, let us characterize the coupling to GSPs for an 
emitter at frequency u> decaying in the vicinity of a free- 
standing 2D graphene sheet. A sketch of the system 
under study can be seen in Fig. [I] (a). The graphene 
sheet is placed in the x — y plane and has a conductiv- 
ity <y{uS), obtained in the random-phase approximation 
[2"2l |2"3"] . This quantity depends on the chemical poten- 
tial of the graphene sheet, fi, the temperature, which we 
consider to be T — 300 K, and the carriers' scattering 
time r, for which we use a value taken from the theo- 
retical predictions [53] such that E T = h/r = 0.1 meV. 
The emitter, modeled in the point dipole approximation, 
is placed at a distance z from the sheet and has a dipole 
moment p. 

When the emitter is close to the graphene sheet, it can 
decay through three different mechanisms: radiation to 
free-space, excitation of GSPs or coupling to absorption 
losses in the graphene sheet. The emitter's total decay 
rate is proportional to the imaginary part of the Green's 
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tensor of the system, G(f, r, lj), the dipole moment p and 
the free-space momentum ko — uj/c: 
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where u p is the unitary vector in the direction of the 
dipole moment, For the sake of simplicity, we take a 
dipole moment perpendicular to the graphene sheet {p = 
pu z ), and then the relevant component of the tensor is 
G zz [17]. This leads to the following expression for the 
Purcell factor, which is the total decay rate normalized 
to the free-space decay rate (r = fco|p| 2 /(3/T?ieo)): 
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where we integrate over the normalized parallel wave vec- 
tor q = fc||//co, q z = \Jl ~ q 2 is the momentum in the di- 
rection perpendicular to the sheet, with lm(q z ) > 0, and 
r p (q) = —aq z / (aq z + 1) is the reflection coefficient of the 
graphene layer for p-polarization, with a = 2tt<j/c be- 
ing the normalized conductivity. The pole of r p (q) gives 
the dispersion relation of the GSPs propagating in the 
graphene sheet, which appear when Im(a) > 0, i.e., be- 
low a critical frequency Hloq ~ 2/i. The contribution of 
GSPs to the total decay rate can be calculated from the 
pole in r p {q): 
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where q p = \1 — or 2 and q v z = — a^ 1 are the normalized 
momentum components of the GSP. 

The inset panel in Fig. [I] (a) shows the Purcell fac- 
tor (solid red line) at 2.4 THz and for /i = 0.2 eV as a 
function of the emitter-graphene distance z normalized 
to the free-space wavelength, Ao = 124^im. The physi- 
cal parameters [i and Ao were chosen, as shown in Ref. 
[T71 to provide a good compromise in the trade-off be- 
tween confinement and propagation length for the GSPs. 
Three different regions can be identified in the inset panel 
in Fig. [I] (a) according to the decay mechanisms: (i) a 
radiative region at large distances (z > Ao/10 for the cho- 
sen parameters) , where the emitter is far enough from the 
graphene sheet and the total decay rate follows r rad /r 
(dotted blue line), which corresponds to the integration 
of the radiative modes in Eq. [2] (0 < q < 1); (ii) a region 
(Ao/10 > z > Ao/100) where the dominant decay chan- 
nel is the coupling to GSPs and the total decay equals 
T GSP /T (green dashed line); and (iii) a lossy region 
when the emitter is very close to the sheet (z < Ao/100). 
Importantly, and as the figure shows, the decay rate of 
the emitter can be enhanced by several orders of mag- 
nitude. Here we are interested in the plasmonic region, 
where the GSP-contribution to the Purcell factor reaches 
values larger than 100 for the parameters we have cho- 
sen. It is interesting to note that similar values of the 



Purcell factor can be obtained for very thin metal films 
[2 5) when the thickness is much smaller than the skin 
depth, which is challenging from the fabrication point of 
view, as opposed to graphene. For higher frequencies or 
smaller /i, larger Purcell factors in the plasmonic region 
can be obtained in graphene: for instance, T/Tq rj 10 3 
at Ao = 64/im for the same chemical potential. 
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FIG. 1: (a) j3 factor for an emitter at 2.4THz as a func- 
tion of the distance to the graphene sheet, z, normalized 
to the free-space wavelength, Ao, for different values of 
the chemical potential (/i = 0.2, 0.1 and 0.05 eV). Inset 
panel: total decay rate (T/To) and decay rates through 
the plasmonic (r GSP /r ) and radiative(r rad /r ) channels, 
(b) Super/subradiance between two emitters mediated by a 
graphene sheet when the two dipoles interact in transmission 
through it. 7 is plotted as a function of the vertical distance 
of the dipoles to the graphene sheet for the three values of 
the chemical potential \x. The red line shows the vacuum 
interaction j vac ■ 

The parameter that accounts for the efficiency of the 
coupling to GSP, the j3 factor, is defined as the ratio of 
the emitter's decay rate through GSP to its total decay 
rate, j3 = T GSP /T. Fig. [I] (a) studies the possibility of 
tuning /3 with the chemical potential. Three values of /i 
are considered: [i = 0.2, 0.1 and 0.05 eV. For each value 
of fi there is a range of z's where (3 is close to 1, which 
corresponds to the region where the decay rate is dom- 
inated by the plasmonic channel (see inset panel). The 
region of high ft can be dynamically tuned with the chem- 
ical potential, which is in turn controlled by means of an 
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electrostatic gating or a chemical doping. In particular, 
when the chemical potential is decreased to 0.1 eV (green 
line) and 0.05 eV (blue line), the GSP appears at larger q- 
vectors, the GSP is more confined to the graphene sheet, 
and the range of distances where /? is high is narrower. 
The capability of tuning plasmonic properties by means 
of a gate potential is the most important advantage of 
graphene compared to thin metal layers. 



III. SUPERRADIANCE IN TWO-DIMENSIONAL 
GRAPHENE SHEETS 

The efficient and tunable coupling of an emitter to the 
SP modes propagating in a graphene sheet can be used 
to modify the interaction between two emitters placed in 
the vicinity of the 2D sheet, similarly to SPs in metal sur- 
faces [3D] . In order to study the GSP-mediated coupling 
between the two emitters, we introduce the normalized 
decay factor, 7, defined as the ratio of the total decay 
rate of a system where the two emitters interact through 
graphene to the decay rate of two uncoupled emitters in 
the presence of the graphene sheet. Thus, it accounts for 
the modification of the collective decay rate due to the 
presence of the second emitter and reads: 
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The decay rate I\j is the contribution to the decay 
rate of a dipole pi placed at fi due to the presence 
of a dipole pj placed at fj and can be written in 
terms of the Green's function that connects r, to f»: 
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Note that 



for i = j we obtain the decay rate in Eq. [T] i.e., T = Tu. 
The value of 7 characterizes two regimes: when 7 > 1 the 
interaction is enhanced due to the presence of graphene 
and the system is superradiant. Correspondingly, when 
7 < 1, there is an inhibition of the dipole-dipole in- 
teraction and the system is subradiant. Interestingly, a 
graphene sheet allows for two different configurations of 
the emitters: interaction in reflection (both emitters at 
the same side of the sheet) or in transmission (emitters 
placed at opposite sides). 

Let us first consider two emitters interacting through 
the graphene sheet placed at opposite sides of the sheet, 
i.e., in the transmission configuration [see the sketch in 
Fig. [T] (b)]. In order to study the behavior with the 
distance to the sheet, we locate the emitters at the same 
z Qzi\ — \z 2 \) and at rj| = 0. We assume the dipole 
moments to be of the same module but anti-parallel, p\ — 
—pi = pu z - The decay rate related to the interaction 
between the emitters, Ti2, is needed to determine 7 and 
is related to the transmission coefficient: 
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FIG. 2: Tuning super- and subradiance: 7 factor as a function 
of the in-plane distance r between two dipoles in the reflection 
configuration for two different frequencies, (a) At v = 2.4 
THz, when the dipoles are placed at Z\ = 2/im (fi — 0.98, 
T/r = 105), z 2 = 15/oti (/3 = 0.55, r/r = 2) and 23 = 
50^m (/3 = 0.01, r/r = 0.92). The free-space wavelength 
is Ao = 124/im, the plasmon wavelength A p = Ao/3.5 and 
the propagation length L p = 14.9A P . The dashed black line 
corresponds to the free-space interaction and the dashed gray 
line to Eq. [F]for high /3. (b) At v = 7.4THz, when the dipoles 
are at z = 0.2/oti, with (3 = 0.98 and Tn/To = 2250. The 
free-space wavelength is Ao = il/im, the plasmon wavelength 
is A p = Ao/10 and the propagation length is L v = 20A P . The 
dashed black line shows the free-space 7 and the gray line the 
decay of the interaction. 



where the — sign comes from the fact that the dipoles 
are anti-parallel and the transmission coefficient is t(q) = 
l/(aq z + 1). The 7 factor as a function of z/X is plot- 
ted in Fig. [I] (b) for the same parameters used in panel 
(a). The red line shows 7 when the emitters are placed 
in free-space and interact only through radiative modes: 
7vac = at z << A Q because the opposite phase of 
the dipole moments inhibits the radiation and when z 
increases 7vac oscillates around 1 with Ao. When the 
graphene sheet is present, the interaction between the 
emitters is strongly modified at the subwavelength scale. 
In the limit of large z, in correspondence with the dis- 
tances where /3 ~ in panel (a), the emitters couple via 
radiative modes and 7 approaches 7vac- O n the other 
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hand, in the range of z where the plasmonic coupling 
between the emitters starts to dominate, ^ 0, 7 devi- 
ates from 7vac, and, as the distance between the emitters 
and the sheet decreases, the system turns from subradi- 
ant to superradiant. For each value of /1 (0.2, 0.1 and 
0.05 eV), the value of z/Xo where /3 starts to grow from 
to 1, is the onset of the separation between 7 and 7vac- 
Thus, the superradiant regime, controlled by high cou- 
pling to GSP, can be tuned by means of fi. In the limit 
z << A , where (3 = again and losses dominate, the 
interaction reaches 7 = 2, in contrast to the free-space 
value in this limit, 7vac = 0. The reason for this is that, 
in this limit, the integrals in Eq. [2] and Eq. [5] are dom- 
inated by the contributions coming from large q, where 
9z = and r n is controlled by — Im[r(q)] and Tf 2 by 
—Im[t(q)]. Since the imaginary part of both coefficients 
is the same, Tn = Tj 2 and thus 7 = 2. It is interesting 
to note that this sign change comes from the continuity 
conditions of the electromagnetic fields. 

Let us now study how the interaction evolves with the 
in-plane distance between the emitters, r = |ri|i — r 2 ||l> 
where rji = (x,y). We take two dipoles interacting 
through the graphene sheet in the reflection configura- 
tion, as sketched in Fig. [2] (a). In this case we place both 
of them at the same distance z to the sheet, separated 
by an in-place distance r, and with dipole moments of 
the same magnitude, parallel and pointing in the vertical 
direction, V\ = Vi = pu z - In order to determine the 7 
factor we need the interaction decay rate in reflection: 
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where the in-plane dependence is given by the zeroth- 
order Bessel function, Jo. The 7 factor as a function of 
r/A is plotted in Fig. [2] (a) at 2.4THz and [i = 0.2 
eV when the dipoles are at three different separations 
to the sheet, z±, Zi and Z3. First, for z\ — 2/im, we 
know from Fig. [I] (a) that the coupling to GSPs is very 
efficient, /3 = 0.98, and Y )Yq = 105. In this situation, the 
interaction between the two emitters is mediated by GSP, 
and, consequently, the length scale of the interaction is 
controlled by A p = A /3.5 as opposed to the free-space 
interaction, dominated by A (dashed line). When /3 ~ 1, 
an analytical expression for 7 can be obtained in the pole 
approximation: 



7 = 1 + l3J (k q P r) 
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where L p is the propagation length of the GSP, given by 
L p = Ao/[27rim(<7p)], and equal to 14.9A P in this case. 
As it can be seen in the plot, the analytical (gray dashed 
line) and exact (red line) calculations of 7 coincide. Since 
the propagation length is large enough, the decay length 
of the interaction is then given by the Bessel function, 
that decays as \f2j[wr). This is a dimensionality factor, 
coming from the fact that the 3D interaction in free-space 



is confined to the 2D graphene sheet. When the distance 
to the graphene sheet is increased, (3 decreases and 7 de- 
viates from the analytical expression. For Z2 — 15/xm 
the /3 factor is 0.55 and the shape of 7 reflects the fact 
that the emitter decays both to GSP and radiatively. Fi- 
nally, when the distance to the sheet is large enough to 
have P « 0, such as z = 50/im, the vacuum interaction is 
recovered. Therefore, our results show that a larger inter- 
action length scale and a modification of the super- and 
subradiant regimes can be achieved in a subwavelength 
scale for the appropriate choice of parameters. 

The tunability of graphene enables us to reach a regime 
where the interaction between the two emitters can be 
controlled at very deep subwavelength scales. Although 
the tuning can also be done via /x, here we show a situa- 
tion where the tuning parameter is the frequency. When 
the two emitters interact in reflection [see Fig. [2] (b)] at 
7.4 THz and /.i = 0.2 eV, 7 (red line) is very different from 
the one corresponding to the free-space situation (black 
dashed line). Increasing the frequency while maintain- 
ing the chemical potential, results in a larger momentum 
for the GSP, q p , which leads to a tighter confinement as 
well as a reduction in the propagation length. Thus, the 
interaction varies in a X p — Ao/10 scale, as opposed to 
the free-space interaction, dominated by Ao [Fig. [2] (b) 
shows both scales r/Ao and r/\ p ]. The decay of the 7 
factor (gray line), is given by the square-root decay char- 
acteristic of 2D interactions for distances shorter than 
the propagation length, that in this case is L p — 1.85Ao 
(19.2A P ). For larger distances, losses start to dominate 
and the interaction decays exponentially, according to 
Eq.[7} 



IV. SUPERRADIANCE IN ONE-DIMENSIONAL 
GRAPHENE RIBBONS 



For completeness, we have also considered the possi- 
bility of confining GSPs in ID graphene ribbons, which 
could provide a platform for long distance entanglement 
between two emitters, as proposed in Ref. [2JJ Compared 
to the GSPs propagating in a 2D graphene sheet at the 
same frequency, the ribbon-GSPs have a higher q vector, 
thus it is more tightly confined to the graphene layer. 
On the other hand, while the coupling of two emitters 
mediated by 2D-GSPs is dominated by the dimensional- 
ity factor, e~ r / Lp /y/r (see Fig. 2), it is expected that in 
the case of lD-GSPs this coupling will decay as e~ x l Lp , 
enabling long range interaction between the emitters pro- 
vided L p is long enough. Let us consider a free-standing 
graphene ribbon of width A at \y\ < A/2, placed at 
z = with its axis along the x direction (see Fig. |3|. 
We take v = 2.4 THz and the fundamental mode of a 
ribbon of width A = 5/im, that originates from the hy- 
bridization of two edge modes and has even parity of 
E z with respect to the ribbon axis [21]. The field pro- 



5 



Re(£-) 



(a) 



graphene ribbon 



(b) 

graphene ribbon 



(c 



A=5um=X /25 
z=3.12um=y40 




0.5 1 



FIG. 3: Interaction mediated by graphene ribbons. (a,b) 
Electric field profile for a dipole decaying to the ribbon-GSP. 
The dipole is placed at x = 0, y — and z — Ao/40 in panel 
(a) and at z — Ao/10 in panel (b) [the same would be ob- 
tained for 2 > Ao/10]. (c) Super-radiance mediated by the 
ribbon-GSP mode shown in panel (a). The green line shows 
the exponential decay of the interaction. 



file of the fundamental mode, obtained by means of the 
Finite Element Method (COMSOL software), is shown 
for two situations, where the emitter is placed at dis- 
tances, z = 3.12/im and z = 12.4/zm, as shown in Fig. [3] 
(a,b). For a distance to the ribbon of 3.12/im, i.e. Ao/40, 
P m 1, the GSP mode is excited with a very high effi- 
ciency and its field structure is clearly seen in panel (a). 
On the other hand, when the emitter is not sufficiently 
close to the ribbon, /3 << 1 and it couples mostly to 
radiation, as can be seen in panel (b) for an emitter at 
z = Ao/10 where the field snapshot virtually coincides 
with an spherical wave. For the case with /3 « 1, Fig. [3] 
(c) shows 7 (blue line) for two emitters interacting in re- 
flection through the ribbon-GSP. Similar to the 2D GSP, 
a subwavelength modification of the interaction can be 
achieved. However, since the ribbon-GSP is much more 
confined (A p = A /6.1), this modification can be achieved 
in a shorter length scale. Moreover, propagation in ID 
allows for a longer interaction range, with an exponen- 



tial decay given by a propagation length L p — 20 X p . Our 
results demonstrate that graphene ribbons could be used 
to control the length scale of the interaction between two 
emitters thanks to the efficient excitation of GSPs. 



When analyzing graphene ribbons it is worth studying 
the transition across the the ribbon, i.e., how the cou- 
pling between two emitters is affected by the presence of 
edges (see Fig. Eb. For this reason, we consider the evo- 
lution of the interaction in the reflection configuration, 
the same as in Fig. [3] (c), and also plotting the evolu- 
tion of j(x) as a function of the lateral distance from 
the center of the ribbon (see Fig. [4] (a)). As sketched 
in the inset panel, wc displace both emitters perpendic- 
ularly to the ribbon axis (along y), from the ribbon's 
center y = and passing through the ribbon's edge at 
y = A/2, up to y = 2A. As the figure shows, 7(1) evolves 
from reflecting a high coupling to GSP at y = 0, to fol- 
lowing the free-space interaction at y = 2A, i.e., /3 goes 
from w 1 to ss in a length scale of the order of 2 A. 
Additionally, we also study the evolution of the interac- 
tion between two emitters with opposite dipole moments 
in the transmission configuration for two situations and 
plot 7 at x — as a function of y in panel (b). First, 
we displace the emitters simultaneously such that both 
of them are at the same y [see sketch (1)]. At y = 0, the 
emitters couple to the ribbon-GSP and the system is in a 
superradiant state (red line), as opposed to the situation 
in free space, that is subradiant (inset panel, red line). 
When the emitters are displaced from the ribbon's cen- 
ter, but are still on top of the ribbon, i.e. \y\ < A/2, 7 
is only slightly modified. Once the dipoles pass the rib- 
bon's edge, the coupling to the ID GSP is reduced, and 
therefore 7 decreases. Subradiance is quickly reached, ap- 
proaching the free-space value, 7 = 0, which is achieved 
when the emitters are placed at a distance 2A from the 
ribbon's center. In the second situation, one emitter is 
kept at the ribbon's center and the other is displaced 
perpendicularly to the ribbon's axis [see sketch (2)]. In 
this case, the system is always superradiant for the dis- 
tances considered: 7 starts at 2 and approaches 1, while 
in free space 7 is of the order of 0.05 at y = A (blue line 
in the inset panel). The reason for this lies on the fact 
that the emitter that is fixed always couples to the GSP. 
Remarkably, with only one emitter efficiently coupled to 
the ribbon-GSP, the interaction between both emitters is 
very different to the vacuum case. Our results for both 
configurations (reflection and transmission) demonstrate 
that the coupling between emitters mediated by 1D-GSP 
is very insensitive to the lateral displacement and that 
the effective lateral extension of these 1D-GSP is of the 
order of A/2 measured from the edge's ribbon. 
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FIG. 4: Dependence of the dipole-dipole coupling in graphene ribbons on the lateral separation from the center of the ribbon, 
y — 0. (a) 7 factor between two dipoles in the reflection configuration as a function of x for several lateral separations from 
the center of the ribbon (blue lines) . The scale in each sub-panel is between and 2. The red line shows the corresponding 
interaction in vacuum. The inset panel shows a top view of the structure, (b) 7 factor as a function of the lateral separation 
y in the transmission configuration for two cases. First (red line), 7 factor when the two dipoles are displaced simultaneously, 
as shown in sketch (1). The dots correspond to numerical simulations and the continuous line is a guide-to-the-eye. Second 
(blue line), 7 factor (from simulations) when the upper dipole is fixed at the center of the ribbon and the lower is displaced, as 
shown in sketch (2). The inset panel shows the corresponding interaction when the dipoles are placed in vacuum. 



V. CONCLUSION 

In conclusion, we have studied the tailoring of the 
interaction between two emitters mediated by surface 
plasmon modes in a graphene sheet. We have shown 
that within a certain range of distances to the graphene 
sheet, the decay rate of one emitter can be fully dom- 
inated by the graphene surface plasmons. Due to this 
efficient coupling, the enhancement of the decay rate 
of the emitter, or Purcell factor, can be enhanced by 
several orders of magnitude. The interaction between 
two emitters mediated by the graphene plasmons in two- 
dimensional graphene sheets can thus be controlled at 
a subwavelength scale and can be tuned by means of 
external parameters. We have studied the appearance 
of the super- and subradiant regimes, both in the reflec- 
tion and transmission configurations. Additionally, when 
the interaction is confined to one-dimension in graphene 
ribbons, a longer interaction range between the emitters 
and a very deep subwavelength control of the interaction 
can be achieved. Here, the lateral confinement leads to 
much higher enhancement factors, very deep subwave- 
length length scales for the coupling and longer interac- 
tion ranges. By considering the lateral displacement of 
the emitters from the ribbon's axis, we have also shown 
that the coupling to the graphene surface plasmons sup- 
ported by the ribbon is very robust. Our results show 
that both graphene sheets and graphene ribbons can be 



used as efficient platforms to modify the interaction be- 
tween two emitters when they are placed in their vicinity. 
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